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ABSTRACT 

We perform global time-dependent simulations of an accretion disc around a young 
stellar object with a dead zone (a region where the magneto-rotational instability 
cannot drive turbulence because the material is not sufhciently ionised). For infall 
accretion rates on to the disc of around 10"'' Mq yr~^, dead zones occur if the critical 
magnetic Reynolds number is larger than about 10*. We model the collapse of a 
molecular gas cloud. At early times when the infall accretion rate is high, the disc 
is thermally ionised and fully turbulent. However, as the infall accretion rate drops, 
a dead zone may form if the critical magnetic Reynolds number is sufficiently large, 
otherwise the disc remains fully turbulent. With a dead zone the disc can become 
unstable to the gravo- magneto instability. The mass of the star grows in large accretion 
outbursts that may explain FU Orionis events. At late times there is not sufficient mass 
in the disc for outbursts to occur but the dead zone becomes even more prominent as 
the disc cools. Large inner dead zones in the later stages of disc evolution may help 
to explain observations of transition discs with an inner hole. 

Key words: accretion, accretion discs - planetary systems: protoplanetary discs - 
stars: pre-main-sequence - circumstellar matter 



1 INTRODUCTION 

The core of a protostellar molecular cloud collapses to form a 
low-mass protostar with a dis c on a timescale of a few 10^ yr 
l|Shu. Adams fc Lizand[l987h . Angular momentum is trans- 
ported outwards through the disc by turbule nce, thus allow - 
ing material to accrete on to the central star (|Pringlelll98ll ) . 
Turbulence in accretion dis cs may be driven by the m agneto- 
rotational instability (MRI. lBalbus fc Hawley||l99ll) but can 
be suppressed by a low ionisation fraction ( Gammij 1 19961 : 
iGammie fc MenoulllQQSl ). The inner parts of the disc may 
be hot enough to be thermally ionised and fully MRI ac- 
tive (turbulent). Further out, where the disc is cooler, a 
dead zone forms at the midplane. In this region the MRI 
cannot drive turbulence. The surface layers of the disc are 
ionised by cosmic r ays or X-rays from the star and may re- 
main turbulent (e.g. lSano et al.l2000l : lMatsumura fc Pudritj 
I2OO3I ). A second form of turbulence may also be driven by 
gravitational ins t ability if the disc bec omes massive enough 
l|Paczvnski|[l97i : iLodato fc Ric3l2004h . 

Previously, time-dependent numerical disc simulations 
have used a simple description for the extent of the dead 
zone llArmitage. Livio fc Pringld [2OO1I : IZhu et al.1 l2010al : 
IZhu. Hartmann fc Gammie 2010bt ). They assume that cos- 
mic rays penetrate a constant surface density in the disc 



surfaces and are attenuated expone ntially with a stopping 
depth of less than 100gcm~^ (e.g . lUmcbavashi fc Nakand 
■ 19811: IZhu. Hartmann fc Gammidf [2009 : Terquem ,20081 : 
iMatsumura. Pudritz fc Thommed |2009| )~ A more realistic 
way to determine the dead zone extent is with a mag- 
netic Reynolds number. Turbulence is suppressed if the 
magnetic Reynol ds number is l ess than some critical value 
Rem < Rem. (|Fleming||2000h. However, as we discussed 
in a preceding paper ( Martin et al.ll201^ . MLLPl), the criti- 
cal value is uncertain. The magnetohydrodynamical (MHD) 
simulations of .Flemin g (2000) suggest it is of the order of 
10* if there is no magnetic flux through the disc, but it may 
be close to 10 if ther e is a m agn etic flux. The linear stabil- 
ity analysis of lWardld (|l999l ) and lBalbus fc Tergueml (|200lD 
suggest it may be as small as 1, but conditions for linear sta- 
bility are n ot necessarily the same as those for maintaining 
turbulence (|Balbus fc Hawlevll200oh . In MLLPl we showed 
that a constant surface density in the active layer may only 
be a good approximation for a disc with a low critical mag- 
netic Reynolds number of around 1. Even then, the surface 
density of the layer may be much larger than the constant 
values that have been assumed, that were less than about 
200 gcm-^ 

The luminosity of the protostar depends upon the ac- 
cretion rate through the disc. The typical bolometric lumi- 
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nosity of protostars is much smaller than would be expected 
from the observed infall rate and duration of the protostel- 
lar phase (|Kenvon et al.|[r990l ). Time-dependent accretion is 
thought to be a solution to this luminosity problem. Ob- 
servational evidence for this scenario has been seen in FU 
Orionis systems that show outbursts with a peak accretion 
rate of a round 10~"^ yr~^ lasting a timescale o f years to 
decades (|Herbielll977l : iHartmann fc Kenvonlll996l '). A plau- 
sible explanation for FU Orionis outbursts is a disc insta- 
bility resulting from the formation of a dead zone, known 
as the gravo-magneto instability. The MRI turbulent ac- 
tive layers of the disc supply material to the dead zone 
that builds up until it becomes self-gravitating. Self-gravity 
drives a second type of turbulence that can dissipate suf- 
ficient heat for the MRI to be triggered. The now very 
high surface density disc flows on to the central star in 
an accretion outburst. The disc then cools, the dead zone 
re-forms and the cycle repeats (lArmitage. Livio fc Pringld 



l200ll : IZhu. Hartmann fc Gam mi j'2009''l . This process can be 
understood as a limit cycle (|Martin fc L ubow 2011a). There 
are two steady state solutions, one with a fully turbulent disc 
and a second that is self-gravitating. There is a range of ac- 
cretion rates for which no steady state disc solution exists. 
The outbursts can be explained as transitions between the 
two steady solutions in a state diagram that plots the ac- 
cretion rate through the disc against the surface density at 
a fixed radius. 

In the later stages of disc evolution, the infall accretion 
slows. There is not enough mass flowing through the disc 
for the gravo-magneto instability to operate on a reasonable 
timescale yet there may still be a dead zone present. Ob- 
servations of protostellar disc s show that about fiv e to ten 
percent are transition discs (|Skrutscki et al.lll990l ). These 
have no significant near-IR excess compared to younger 
discs and are thought to have a hole in the inner regions 
up to sizes of seve ral tens of AU (e.g. ICalvert et al.l [20021 : 
iForrest et J] 120041 '). There is still some gas accretion on to 
the star but no dust is observed there. There have been 
sev eral models sugge sted to explain these discs such as plan- 
ets ([Rice et al.llioOa '). disc evolution JPullemond fc DominikI 
2005!'). disc clearing ((ciarke. Gendrin fc SotomavoJ I2OOII : 
Alexander. Clarke fc Pringlell200^ and the presence of com- 
pani ons ('Je nsen fc Mathieul 1199*^ ). [Alexander fc Armitagel 
(|2009) suggest that young discs contain a planet whereas 
the older ones may be undergoing disc clearing. We consider 
an additional possibility that transition discs may contain 
a large dead zone that causes a sharp transition in the sur- 
face density of the disc and allows only a small amount of 
accretion on to the central star. 

In Section [2] we describe the layered disc model. In Sec- 
tion [3] we find fully turbulent steady state disc solutions 
with a constant infall accretion rate and consider parame- 
ters for which a dead zone can form. In Section U we solve 
the time-dependent layered accretion disc equations and find 
the evolution of a system with a dead zone and a constant 
infall accretion rate. In Section [S] we model the infall ac- 
cretion rate from a collapsing gas cloud and consider disc 
parameters necessary for the gravo-magneto instability to 
explain observed FU Orionis outbursts. In Section[S]we dis- 
cuss implications of a dead zone model for transition disc 
observations. 



2 LAYERED DISC MODEL 

In an accretion disc, material orbits the central mass, M, 
at Keplerian v elocity at radius R with angular velocit y 
n = -JGMJW (|Lvnden-Bell fc PringlellT97i : IPringlelll98lh . 
We use a one-dimensional layered disc model to follow the 
time evolution of the total surface density, E, and mid- 
plan e temperature, Tc as described by equations 1-8 in 
Mar tin fc Lubowl (|2011al ) but we do not assume a constant 
surface density in the active layer. Instead, we use the an- 
alytic approximations for the active layer surface density 
given in equations 26 and 27 in MLLPl. These are approx- 
imations to the solution of the equation Rem = i?eM,crit, 
assuming that cosmic rays are the dominant source of ionisa- 
tion. We note that this description assumes Ohmic resistiv- 
ity is the dominant non-ideal MHD effect. However, recent 
works suggest that ambipolar diffu sion and the Hall effect 
may d etermine the active la yer (e.g. Perez-Becker fc Chiand 
l2011a .b: Bai fc Stone 20ll ) . IWardle fc SalmeronI (|201ll ) find 
the Ohmic resistivity term provides an average value for 
the active layer surface density for a range of vertical mag- 
netic fields. Because Ohmic resistivity also has th e advan- 
tage of being independent of the magnetic field (|Fleming| 
l2000h we consider only this effect here. We note that other 
sources of turbulence, such as the baroclinic inst ability (e.g. 
i Klahr fc Bodenheimeijl2003l : lLvra fc Klahij|201 j ). have been 
suggested to generate viscosity within the dead zone. How- 
ever, in this work we we assume there is no turbulence within 
the dead zone. 

In MLLPl we found that the dead zone extent is not sig- 
nificantly affected by variation in the metallicity and hence 
we use the high metallicity limit for the electron fraction (see 
equation 19 in MLLPl) and consider changes in the criti- 
cal magnetic Reynolds number. The active layer has a sur- 
face density Em, temperature Tm and viscosity parametrised 
with the a prescription by 

y = aCmH (1) 

feha kura fc Sunvaevlll973l ). where the sound speed in the 
active layer is 



(2) 



where k is the Boltzmann constant and /irrih is the mean 
molecular mass. We approximate the disc scale height by 
H — Cm/ii and so 



(3) 



The lShakura fc Sunvae Q parameter still has some 

uncertainty. Q Numerical MHD simulations find a — 0.01 



^ We note that it may be more realistic to take u = acmh where 
h is the local scale height in the layer. The scale height in the 
active layer should he h = min{H, p/^) at the base of the active 
layer and so typically h <^ H. As mass builds up in the dead 
zone, the magnetic surface layer will lie further above the disc 
midplane, resulting in higher local gravity and a smaller scale 
height. This means that the a parameter we take is effectively 
increased to a^g = {H/h)a. For example if we have a = 0.01 
and H/h = 10 then the effective parameter in the active layer is 
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Figure 1. Steady state disc solutions that are fully turbulent with a = 0.1 (left) and a = 0.01 (right). The solid lines show the total 
surface density. The upper lines in each plot show a disc with an infall accretion rate of M^^f^n = lO~^M0yr~^ the middle lines 
A^infall = lO^'^MQyr"^ and the lower lines Minfaii = lO~*M0yr~^. The dashed lines show the surface density that is turbulent for 
Rbm crit = 10* a-nd the dotted lines for -Rem ^rit = W^- The dotted lines are not visible in all cases because the disc is fully turbulent 
and the active surface density is equal to the total surface density. However, if a dead zone exists in the disc, then the disc can not be 
in steady state. 



l|Brandenburg et al.l [l995l : IStone et al.1 Il996l) but observa- 
tions suggest a ~ 0.1 - 0.4 (jKing. Pringle fc Livid lioOTl ) . 
If the disc is fully turbulent the single layer has Em = E 
and Tm = Tc. However, where a dead zone exists the disc 
has two layers and the complementary region is defined with 
surface density Eg = E — Em and temperature Tc. The dead 
layer may become self-gravitating if Q < Qcrit, where the 
Toomre parameter is 



^ ttGE' 

the sound speed at the midplane is 



(4) 



(5) 



and we take Qcrit ~ 2. Self-gravity driv es a second effec- 
tive vis cosity in the disc as explained in iMartin fc Lubowl 
(|2011al ). 

In MLLPl, we found that the active layer surface den- 
sity may be significantly smaller than that assumed by the 
constant prescription, especially at the inner edge of the 
dead zone. The active layer may become optically thin and 
this affects the relationship between the layer temperatures. 
The optical depth of the turbulent layer is 



3 ^rp \ ^ni 

"^m — g^V-tmiPmJ ^ ' 



(6) 



where k is the opacity and the density is approximated as 



larger, a^g = 0.1. This should be investigated further in future 
work. 



Pm = Em/2-ffm with scale height Hm ~ Cm/f2. The optical 
depth within the complementary region, if it exists, is 



3 IT 



Pcj — , 



(7) 



where the density is pc = Eg/2_ffg and the scale height is 
Hg = Cg/Sl. The total optical depth of the disc is 

r = rm + rg, (8) 

the sum of the optical depths of the two layers. 

In this work we consider two different opacity prescrip- 
tions. The first is a simple power law in temperature 



K(r) = 0.02T''®cm^ 



(9) 



(jArmitaee. Livio fc PringlellioOl] : iMartin fc Lubowl 1201 lal l. 

This is useful for finding analytic solutions because it does 
not depe nd on the density o f the disc. It is the opacity in the 
tables of iBell fc LinI (|l994l ) around the t emperature where 
the MRI is trig gered of around 800 K (see lUmebayashilll983l : 
iGammie 19961). The second is the f u ll opa city shown in ta- 



Zhu. Hartmann fc Gammid (|2009l ) that depends on 



ble 1 of 

both the temperature and the density of the disc. 

If the active layer is optically thick, the opacity 
the temperatures in the layers are found with 

m — ' m c ; 

where Te is the surface temperature, and 

lEm^l Tcr . 



> 1, 



(10) 



o 



However, if the active layer is optically thin, 
Tm = To and 



(11) 
< 1, we take 

(12) 
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The active layer may be optically thin in the inner parts of 
the disc where the layer is thinnest. 

With a large active layer surface density the inner parts 
of the disc remain hot enough for thermal ionisation. How- 
ever, when the dead zone is determined by a critical mag- 
netic Reynolds number the active layer can become very thin 
and the inner parts of the disc provide little accretional heat- 
ing. In this case, heating from the star may dominate the ac- 
cretional heating. We consider this further in Section r4.3l but 
find it makes little difference except for radii R < 0.1 AU. 



3 ANALYTICAL STEADY STATE DISCS 

For sufficiently small critical magnetic Reynolds numbers a 
steady state disc solution exists that is fully turbulent, thus 
Em = E everywhere. There is no complementary region. 
Eg — 0, and the temperatures are related with Tn, 

n4 



Tc and 

r," = rT^ (13) 
For a fully turbulent disc, the surface density is found with 



3^ 



1 - 



-Rin 



(14) 



l|Pringle|[l98ll ). where M = Minfaii is the steady accretion 
rate through the disc. There is a zero torque inner boundary 
condition applied at R = i?in = SR© where the mass falls 
freely on to the star. Similarly, the surface temperature is 
found with 



8-K 



Rin 



(15) 



l|Pringldll98lh . In this section, we use the simplified power 
law opacity given in equation Q so that we can easily solve 
equations (|13|l . (|14|) and (|15p analytically for the surface 
density and temperature in the disc. With this fully turbu- 
lent disc solution, we check if the solution has a dead zone 
(we solve Rcm = iieM.crit). If it does, then the steady state 
does not exist and we expect time dependent accretion (see 
Section IH). 

In Fig. [T] we show the steady state surface density of 
the disc around a star of mass M = 1 Mq for varying in- 
fall accretion rates with critical magnetic Reynolds number 
-ReM.crit = 10^ and 10* with a = 0.1 on the left and a = 0.01 
on the right. For larger a, the steady state surface density 
is smaller. Unless the accretion rate is very small, there is 
no dead zone for -ReM.crit ^ 10*. 

In Fig. [2] we show the corresponding disc temperatures, 
for the fully turbulent disc, that decrease with the infall 
accretion rate. Because the active layer surface density de- 
creases with temperature (MLLPl) the active layer surface 
density also decreases with the infall accretion rate. In a 
time dependent disc, this steady solution can only exist if 
the disc is fully turbulent. 

In MLLPl we found that the active layer surface den- 
sity is only sensitive to the temperature and not the total 
surface density. The temperature in the steady state disc 
found here is higher than often consi d ered in static discs (e.g. 
Fromang. Terguem fc Balbui l2002l : iMatsumura fc Pudrit j 
20031 ). The new result is that for time evolving discs, dead 
zones exist only for large critical magnetic Reynolds num- 
ber. In order to investigate discs with large critical magnetic 
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Figure 2. The temperature of the steady state discs shown 
in Fig. [1] with a = 0.1 for Minf^ii = 10-''MQyr-i (solid 
line), Minfaii = lO"'^ Mq yr~^ (short-dashed line) and Minfall = 
10~^ Mq yr~^ (long dashed line). The dotted line shows a = 0.01 
and Minfall = 10-^ Mq yr-\ 



Reynolds numbers, and hence dead zones, we must solve the 
time-dependent accretion disc equations. 



4 TIME-DEPENDENT SIMULATIONS 

We solve numerically the time-dependent layered accre- 
tion disc equations with a one-dimensional code similar to 
iMartin fc Lubowl ()2011al ). The disc around a star of mass 
M = 1 Mq extends from Rj^ = 5 Rq up to -Rout = 40 AU 
(see e .g. lArmitage. Livio fc Pringlj l200ll : iMartin fc Lubowl 
l2011al lbh. The grid has 200 points distributed uniformly in 
log R so that at each radius the grid point separation is 
AR/R = 0.037. Material is continuously added to the disc 
at a radius of -Radd = 35 AU at a rate Minfaii . The inner edge 
of the disc has a zero torque inner boundary condition and 
the outer edge has a zero radial velocity boundary condition 
that prevents outward flow. The initial conditions do not 
have any effect on the outcome because we evolve the disc 
until it reaches a steady state or a limit cycle in the accretion 
rate. We consider both the simple power law opacity given 
in equation ([9|l and also a full op acity given in table 1 of 
IZhu. Hartmami fc Gammid (|2009| ) in our numerical models. 



4.1 Small Critical Magnetic Reynolds Number 

Only with large critical magnetic Reynolds numbers does a 
dead zone exist in the steady state disc. However, our ana- 
lytical models in the previous section used only the simple 
power law opacity. In this Section, we numerically evolve the 
disc until it reaches a steady state (with no dead zone) in 
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R/AU R/AU 

Figure 3. The fully turbulent steady state disc surface density (left) and temperatures (right) with an infall accretion rate Afinfaii = 
10~^ M0 yr~^ and a = 0.01. The solid lines show the analytical model given in Section |3] with the simple opacity law. The dotted lines 
show the corresponding numerical model with the simple opacity law (this is almost identical to the analytical solution). The dashed 
lines show the numerical model with the full opacity table. For the temperature plot we show the midplane temperature, Tc, (upper 
three lines) and surface temperature, Te (lower three lines). The surface temperatures are almost identical for all three models. 



order to examine properties of a disc with the full opacity 
prescription. 

In Fig. [3] we show the surface density of the fully tur- 
bulent steady state disc with an infall accretion rate of 
Minfaii = 10~^ Mq yr~^ and a = 0.01. The simple power law 
opacity approximates well the analytical solutions shown. In 
the outer parts of the disc, the numerical solutions have a 
zero radial velocity condition imposed that alters the solu- 
tion slightly from the analytical one. We also show a simu- 
lation for a disc with the full opacity table. The inner parts 
of the disc are significantly more massive than with the sim- 
ple opacity law. However, these parts are hot enough to be 
thermally ionised. 

In Fig. [3] we also show the midplane and surface tem- 
peratures of the steady state disc. The numerical simulation 
with the simple opacity law approximates well the analytical 
steady state solution. The simulation with the full opacity 
table shows that the inner parts of the disc may be slightly 
cooler compared with the simple opacity law. However, fur- 
ther out in the disc where dead zone formation occurs the 
power law is a reasonable approximation to the full opac- 
ity. Hence, even though we have taken a simple opacity law 
in the previous section, it does not change our conclusion 
significantly that dead zones will only form for large criti- 
cal magnetic Reynolds numbers iieM.crit ^ 10* unless the 
accretion rate is very small (we have also verified this nu- 
merically) . 



4.2 Large Critical Magnetic Reynolds Number 

For sufficiently large critical magnetic Reynolds number, 
the steady state disc solution contains a dead zone. Hence, 
it cannot be in steady state and we expect the disc 



to be unstable to the gravo-magneto instability. This re- 
sults in ou tbursts in the accretion rate on to the cen- 
tral object (IGammi 3ll996l:lArmitage. Livio fc Pringle|[200ll : 



IZhu. Hartmann fc Gammiel 20091 ). We use the analvtic ap- 
proximation for the active layer surface density given in 
equations 26 and 27 in MLLPl that is valid for 7?eM,crit 
100. We choose -ReM.crit = 10**, Minfaii = 10~^ Mq yr"* and 
a = 0.01 so that we can compare with previous work with a 
constant surface density in the active layer. We use the sim- 
ple power law opacity because we have seen in Section 14.11 
that it makes little difference except to the hottest parts of 
the disc which are already thermally ionised and fully tur- 
bulent. 

In Fig. |4] we show the accretion rate on to the central 
star. When the dead zone is present, the accretion rate on 
to the star is extremely low. The accretion outbursts occur 
on a timescale of 1.1 x 10^ yr and a mass of around 0.1 Mq is 
accreted during each. After a few outburs ts the disc reaches 
a limit cycle that repeats itsel f . Follo wing iMartin fc Lubowl 
lioili) and lLubow fc Martini (|2012| ). in Fig. [5] we show the 
limit cycle in the E — Ms plane at a radius of i? = 1 AU. 
The steady accretion rate Ms corresponds to the steady sur- 
face temperature in equation (|15p . The disc moves in an 
anti-clockwise direction around the cycle. In the dead zone 
branch, the active layer deposits material on to the dead 
zone and as a result, the surface density increases along the 
branch. The accretion rate decreases with increasing surface 
density because the active layer decreases. The disc becomes 
self-gravitating along the gravo-magneto (GM) branch. The 
temperature, and hence accretion rate, begins to rise as the 
surface density continues to increase. The disc eventually 
becomes hot enough for thermal ionisation to take place at 
some radius and the MRI is triggered at the end of the GM 
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Figure 4. The accretion rate on to tlie solar mass central star 
for a disc with Minfaii = 10~^ yr~^, a = 0.01 and Rcm ^rit = 
10*. A mass of 0.1 Mq is accreted on to the central star during 
each outburst. 



branch. The disc moves quickly to the upper MRI branch. 
Now the accretion rate on to the star is much larger than the 
accretion on to the disc and so the accretion rate declines as 
the surface density decreases. The dead zone re-forms, the 
disc quickly transitions down to the dead zone branch and 
the cycle repeats. 

In Fig. [6] we show the surface density and midplane 
temperature of the disc just before an outburst and during 
an outburst. There is very little mass in the inner parts 
of the accretion disc. This is different to the disc with a 
constant (and large) surface density in the active layer. The 
active layer here is very thin and little mass flows around the 
dead zone. Material in the dead zone builds up and peaks in 
the disc where it becomes self-gravitating. The temperatures 
also peak at this radius. The MRI is first triggered at the 
peak and the outburst begins. We have found that the gravo- 
magneto instability mechanism is the same no matter how 
the extent of the dead zone is determined. 



4.3 Irradiation from the star 

When the active layer is very thin, in the innermost parts 
of the disc the accretional heating may be dominated by the 
heating from the star. The flux of radiation is 



where 



air 



0.005 



4 Oirr 



AV J 



-Km)' 



(16) 



(17) 



l|Chiang fc Goldreich|[T997b and we take the radius of the 
star to be 7?* = 5Kq and the temperature T-* — 4000 K. 



Figure 5. The limit cycle in the S — AIs or S — To plane at 
R = I AU for the disc in Fig. |4] The disc moves in an anti- 
clockwise direction around the cycle. The steady accretion rate 
Ms corresponds to the steady state surface temperature given in 
equation l|15|l . 



If the midplane temperature of the disc drops below the 
irradiation temperature 



Tir 



a 



(18) 



then we assume that the disc is isothermal and set Tc = 

When we include this term the inner parts of the disc 
are warmer and as a result, up to a radius of around 0.1 AU 
the disc is thermally ionised. However, because the inner- 
most part of the dead zone still has a very small active layer, 
the accretion rate on to the star remains small. For example, 
for the model in the previous section, without the irradia- 
tion term we find that the accretion rate between outbursts 
is around 10~^^ Mq yr~^ whereas with the effects of irradi- 
ation this increases to around 10~^^MQyr~^. This is still 
significantly too small to account for the T Tauri accretion 
rates that ar e thought to be around 10~^ Mq yr~^ at an 
age o f 1 Myr (|Valenti. Basri fc Johns|[l993l : iHartmann et al.l 
'iQOSf). For the rest of this work we ignore the heating from 
the star because is only affects a very small part of the disc. 



5 FU ORIONIS OUTBURSTS 

In a collapsing molecular cloud the accretion rate on 
to the disc will not be constant in time. We follow 
lArmitage. Livio fc Pringld l|200ll l who use a simple model 
for the the protostellar accretion history. At early times the 
infall rate is expected to be Mi,jfaii_= Cg/G where Cs is the 
sound speed in the cloud ()Shuiil97711 . For a temperature of 
10 K, this is around Afinfaii ~ lO"'^ Mq yr-\ We assume that 
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Figure 6. Numerical simulations with M = 1 M0, A/infaii = lO"*^ Moyr-\ Rejvi.crit = 10'* and a = 0.01. The surface density (left) and 
midplane temperature (right) in the disc just before an outburst (solid lines and shaded region), at the start of an an outburst (dotted 
lines) and towards the end of an outburst (dashed lines). Just before the outburst, the dead zone is shown in the shaded region. The 
darker shaded part (labeled GM) shows where the dead zone is self-gravitating. The upper dead zone layer boundary is too close to the 
total surface density to be seen to be different on this scale (the active layer is very thin). During the outburst the whole disc is turbulent. 



the accretion rate declines exponentially in time 

Minfaii = Mexp ^-^^ , (19) 

where we take the initial accretion rate to be Mi = 2 x 
10~^ M0 yr~^. At this time, we begin with an initially fully 
turbulent steady state accretion disc. The free fall timescale 
is 

cloud / 

where pcioud is the cloud densi ty. We take tg — 10* yr 
ijArmitage. Livio &: Prin"gi3l200ll ). 

W ith a star form ation rate predicted bylMiller fc Scalol 
l|l979h . theoretically, iHartmann fc KenvonI (|l985h predict 
that each FU Orionis system should show around ten out- 
bursts during its lifetime. In Fig. [7] we show the accretion 
rate on to the central star as a function of time for the 
infall accretion rate on to the disc given in equation (|19p 
with a — 0.01 and three different critical magnetic Reynolds 
numbers. In all cases, initially the disc is fully turbulent and 
the accretion on to the star proceeds at a rate close to Afinfaii. 
However, once a dead zone forms, the accretion rate drops. 
For a solar mass star with a disc with a = 0.01, the critical 
magnetic Reynolds number must be of the order of a few 
times 10* for a significant number of outbursts to occur. For 
the lowest critical magnetic Reynolds number shown of 10'', 
a dead zone forms once the infall accretion rate has dropped 
to about 3 X 10~^ Mq yr~^. However, at this time there is 
not enough material in the disc for outbursts to occur. 

In Fig. [5] we show the accretion rate on to the star for 
a = 0.1 and -ReM.crit = 5 x 10''. For a higher a we need an 
even larger critical Reynolds number for outbursts to occur. 



We note that the evolution is also dependent on the mass 
of the star. The active layer surface density decreases as the 
mass of the central star increases (MLLPl). Increasing the 
mass of the central star has the same effect as increasing 
the critical magnetic Reynolds number. So, for larger mass 
stars, dead zones will form for a smaller critical magnetic 
Reynolds number. 

If the observed FU Orionis outbursts are driven by 
the gravo-magneto instability, then the critical magnetic 
Reynolds number must be very large, of order a few 10*. 
The outbursts are similar to those previously found with 
a model with constant surface density in the active layer. 
However, we predict a very small accretion rate in between 
outbursts. We should try to observe FU Orionis systems in 
quiescence in order to draw conslusions about the surface 
density in the active layer. It is possible that effects such as 
ambipolar diffusion or the Hall effect could play a role in in- 
creasing the active layer surface density in the inner regions 
and this should be investigated in future work. 



6 TRANSITION DISCS 

At later times of disc evolution, the infall accretion rate on to 
the disc slows and the mass and temperature of the disc both 
decrease. A dead zone determined by a critical magnetic 
Reynolds number becomes more prominent as the disc cools 
(e.g. see Fig. 1). However, with the constant active layer 
surface density previously assumed this is not the case. The 
constant active layer model predicts that when the surface 
density drops below the critical value (of around 200 g cm"'^) 
the whole disc becomes turbulent. The increasingly large 
dead zone found with a critical magnetic Reynolds number 
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Figure 7. The evolution of the accretion rate on to a solar mass 
star. The dashed lines show the exponentially decreasing infall 
accretion rate given in equation I I19II . The solid lines show the 
accretion rate on to the star with a = 0.01, Re^ crit = 10'' (top), 
10"* (middle) and 5 X 10"* (bottom). Initially the infall accretion 
rate is high enough that there is no dead zone. The accretion rate 
on to the star proceeds at a rate similar to that on to the disc. As 
the infall accretion rate drops, the temperature decreases and a 
dead zone forms. Then the accretion proceeds in large outbursts 
with a very small rate in between. At late times there is not 
sufficient material left in the disc for outbursts to occur. When the 
critical magnetic Reynolds number is 10'^, there is not sufficient 
mass left in the disc for outbursts by the time the dead zone has 
formed. 
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Figure 8. As Fig. [7] but with a = 0.1 and ReM,i 



: 5 x 10-* 



has implications for low mass discs in the later stages of 
evolution. 

After the infall on to the disc ends, there is not suffi- 
cient accretion through the disc for it to build up enough 
material for accretion outbursts (see for example top pan- 
nel in Fig. [7|). We have investigated this numerically and 
find the surface density and temperature distributions are 
similar to the pre-outburst distribution shown in Fig. [6] The 
inner parts of the disc, that are mostly dead, are depleted by 
several orders of magnitude compared with the outer parts 
that are self-gravitating. Similarly the temperature is sev- 
eral orders of magnitude smaller in the inner regions. Only 
a small amount of material can flow through the thin active 
layer on top of the dead zone. The depleted inner parts of 



the disc could appear to form a hole at the inside on the 
viscous timescale 
If 

V 



1.41 X 10' 
R 



{—V 

VO.Ol/ 



( T 



VlOOK 



M 



X 



lAU 



(21) 



This is very short compared with the lifetime of the disc. 

A large dead zone could help to explain observed transi- 
tion objects without the need for a photoevaporation model. 
The surface density distribu tion is qualitatively similar to 
that s hown in Figure 1 in IClarke. Gendrin fc Sotomavoij 
(20o3) for their ultra violet switch photoevaporation model. 
However, that model takes a few 10® years to form the in- 
ner hole whereas the dead zone model forms a hole much 
more quickly. We have not included dust grains in our 
models. These suppress the ionisat ion of the disc further 
leadi ng to a larger dead zone (e.g. iBai fc G oodmanl l2009l : 
Perez -Becker fc Chiandl2011al : lBai fc Stonel l2011'). The dead 
zone model should be investigated further in future work 
with respect to observations of transtion discs. 



7 CONCLUSIONS 

We have presented global time-dependent calculations of cir- 
cumstellar discs with a dead zone determined by a critical 
magnetic Reynolds number. For infall accretion rates around 
lO~^M0yr~^, if the critical magnetic Reynolds number is 
larger than about 10^, the disc around a young stellar object 
forms a dead zone. If the critical magnetic Reynolds number 
is smaller, then the disc is fully turbulent. 

If the infall accretion rate is sufficiently high when a 
dead zone forms, the disc becomes unstable to the gravo- 
magneto instability. The outburst behaviour is similar to 
that found previously by assuming a constant active layer 
surface density. If the gravo-magneto instability is respon- 
sible for FU Orionis outbursts, then the critical magnetic 
Reynolds number must be greater than a few 10^. 

A disc model with a dead zone determined by a crit- 
ical magnetic Reynolds number predicts a very low accre- 
tion rate between outbursts that was not predicted with 
the constant active layer surface density assumption. This 
low quiescent rate should be observed both before and after 
FU Orionis outbursts take place. Alternatively, the accre- 
tion rate the model predicts may be increased if effects such 
as ambipolar diffusion or the Hall effect are included. 

The size of a dead zone in a disc increases as the disc 
cools. We suggest that large inner dead zones may help to 
explain observations of the depleted inner regions of transi- 
tion discs. However, there are some limitations to the model 
used here. Some of thes e have already been discussed in 
iMartin fc Lubowl (|2011al ') and MLLPl. There is much more 
work to be done in future to establish criteria for the for- 
mation of dead zones. This is not straightforward because 
it depends on assumptions about the magnetic field and the 
still unresolved problem of the viscosity a parameter. Once 
the criteria are established, the evolution of a protostellar 
disc should be computed, from times of high infall rates to 
late time disc dissipation. 
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